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We present a theoretical model for C1LV3O7: the 1/4- 
depleted square spin- 1/2 Heisenberg model which includes 
both the nearest-neighbor coupling (J) and the next-nearest- 
neighbor coupling (J'), where J and J' are antiferromag- 
netic. Recent experiments of the neutron diffraction by Ha- 
rashina et.al. report the magnetic ordering at low temper- 
atures, which may be called as a stripe phase. It is shown 
that the observed spin structure is not stable in the classi- 
cal theory. By employing the modified spin wave theory, we 
show that the stripe phase is stabilized by the quantum fluc- 
tuations for J'/J > 0.69. In CaVa07, the coupling constants 
are estimated as J ~ J' by comparing the theoretical and 
experimental results. 
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There has been great theoretical and experimental in- 
terest in the physics of quantum phase transition in 
two-dimensional spin systems. Recently, the family 
of calcium- vanadium-oxide, CaV n 02 n +i (n — 2,3,4), 
is added in the list [jlj. These compounds have lay- 
ered structures, where each layer is a regularly depleted 
square lattice. Because V 4+ -ion possesses one d-electron, 
CaV n 02„+i is represented by l/(n + l)-depleted square 
spin-1/2 Heisenberg model. Experimentally, magnetic 
properties of CaV n 02 n +i compounds are quite different 
from one another. A spin gap behavior is observed for 
n = 4, but not for n = 3 [0 |3|. Peculiar geometry of these 
systems together with effects of frustration introduced 
by the next-nearest neighbor coupling result in unusual 
ground state properties. Recently, intensive theoretical 
studies were done for the ground state of CaV^Og [p|-^0[. 
It is now established that some frustration introduced by 
the next-nearest-neighbor coupling is indispensable for 
the realization of the new type of spin disordered state 
(plaquette-RVB state) f§-§. 

In this letter, we study the ground state spin struc- 
ture of CaVaOy. According to the recent experiment 
by Harashina et. al. g], this system shows a magnetic 
ordering at low temperatures below Tn = 23 K. The lat- 
tice structure (1/4-depleted square lattice) of C&V3O7 
together with six-spin unit cell is shown in Fig. 1. The 



nearest- and next-nearest-neighbor coupling constants 
are expressed by J and J', respectively. We assume 
J, J' > because both of them are expected to originate 
from the superexchange through the p-orbital of the oxy- 
gen Q. In this letter, we will show that the structure 
observed experimentally (the stripe phase, Fig. 2(b)) is 
not a stable state but constitutes a saddle-point in the 
classical theory. It is the quantum fluctuations which 
stabilize the stripe phase in the 1 / 4-depleted square lat- 
tice. In many cases the continuous degeneracy existed in 
the classical solutions is lifted by the quantum fluctua- 
tions, order from disorder phenomena H]. However new 
feature here is that a classically unstable state is selected. 

The Hamiltonian of the 1/4-depleted square Heisen- 
berg model is written as 



H = J ^ Sj 

(id) 



n.n.n 



J' \ Si -Sj (J,J' >0), (1) 



(id) 



where Si represents the S = 1/2 spin operator located 
at i-th site. First, we consider the classical energy of the 
Hamiltonian (|l|). In Fig. 2, we illustrate some typical 
spin configurations. Within the six-spin unit cell, the 
classical ground state has the Neel ordering for J' /J < 
1/4 and is changed to the fan structure for J'/J > 1/4 
|§. In the fan phase, 6 = [w - cos -1 ( J/4 J')] /2 and 
<p = W. In the classical theory the transition between 
these states is of the second order. Here, we stress that 
the Neel phase for J'/J > 1/4 and the stripe phase for 
any J'/J are not stable, i.e., saddle-point unstable states. 

In the following, we will consider the effects of quan- 
tum fluctuations by using both the linearized spin wave 
(LSW) theory jH| and the modified spin wave (MSW) 
theory |Q. The MSW theory is applicable to the 
collinear phases, the stripe phase and the Neel phase. 
We express S — 1/2 spin operators by boson operators 
Oj using the Dyson-Maleev transformation, 



S* = I 
2 



ajai, 



Si =4 
i 



"2 + fl . 8 " 
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Sf = (1 - a\ai)a,i, 

for f —spin sites (2) 

Sf = -aj(l - a\a t ), 

= —ai for | —spin sites. (3) 

Then, the spin Hamiltonian (Q) is expressed by the boson 
operators up to quartic terms. We express the sites which 
has three (two) nearest neighbors as A(B)-sites. Clearly, 
the spin reduction due to the quantum fluctuations for A- 
sites and that for B-sites are in general different. Thus, it 



1 



is necessary to introduce two kinds of chemical potentials, 
fi a for bosons at A-sites and fib for bosons at B-sites. Now 
the Hamiltonian (|l|) is represented as 

H = \ H -''J ("-"' • "]"■• "X "'"•) 

(ij)eAF 

+ \ Ji J {- a i a i ~ a j a 3 + a l a 3 + a l a i) 

+ J ij ■ \ a \ ( a ] - a i) a j 

(ij)eAF 

<i,j>eF(T,T) 

+ y^Mi (a] a * _ 2 ) + Classical, ( 4 ) 

where AF (F) represents the summation over the clas- 
sically antiferromagnetic (ferromagnetic) pairs. Jij = J 
(</') for of a (next-)nearest-neighbor pair. It is also 
understood that /i, = fi a i^b) for i £ A-sites (B-sites). 
As the next step, we reduce Eq. (||) to the quadratic 
form by means of a mean- field treatment. For the stripe 
or the Neel phase, three mean fields are introduced. 

By diagonalizing the mean-field Hamiltonian by the 
Bogoliubov transformation in the momentum space, we 
get the magnon-spectra, w£, where k is two-dimensional 
momentum in the magnetic Brillouin zone and a repre- 
sents a branch of the spectra, w£ < (a = 1, 2, • • ■). 
The index a runs from 1 to 3 for the stripe or the Neel 
phase, and from 1 to 6 for the fan phase, reflecting the 
difference on symmetry. Needless to say, they should be 
real and positive for any k and a. At zero temperature, 
the solution of the MSW theory should satisfy the fol- 
lowing relations ; 

<ZI = A, (5) 
1/2 = (4ai) QF + (otoi)BC (* € A, B), (6) 

where A is the non-negative constant. In Eq. (|J), (• • -)qf 
represents the expectation value at zero temperature due 
to zero-point fluctuations, and (■ • -)bc comes from the 
Bose condensation of (k = 0, a = l)-mode, which occurs 
only in case A = 0. The sublattice magnetization for 
A(B)-site, Mmb)i which is derived from the asymptotic 
form of the spin-spin correlation function, is given by 

M A ( B ) = {ala t ) B c for i £ A(B). (7) 

In an ordered phase, there is a Goldstone mode by sym- 
metry requirement and thus the corresponding Bose con- 
densation occurs. In this case, the condition 

A = 0, (oJai)BC > (8) 



should be satisfied. On the other hand, a disordered 
phase solution has a finite excitation gap A for k = 0, 
and no Bose condensation exists. Thus, the condition 

A > 0, (ota i ) B G = (9) 

is satisfied in the disordered phase. In the MSW theory 
we solve self-consistent equations for /i a , fi and the three 
mean- fields with the conditions (||) and (|^) ] (or (^|)). 
Clearly, at the second-order critical point the conditions 
A = and (aja^ec = are satisfied. 

Let us start from the case J = 0, where the prob- 
lem reduces to the one-dimensional chain shown in Fig. 
3. According to Lieb's theorem p5| , the total mag- 
netization of the ground state is rigorously given by 
|(S* ot )| = S ■ |# of A-sites - # of B-sites| = 5A^/3, i.e., 
the ferrimagnetic phase. Here N, which is the total spin 
number of the chain, is assumed to be a multiple of 3. 
Table I shows the ground state properties calculated by 
the LSW, MSW theories and the exact diagonalization 
combined with finite size scaling up to 24 spins (ED). 
Note that these three results satisfy the Lieb's theorem. 
The results by the ED are almost rigorous and relative 
errors are estimated to be less than 10 -3 . It is remark- 
able that even in the system where two kinds of sublattice 
magnetizations (Ma ^ Mb) coexist, results by the MSW 
theory are reliable and are much better than those by the 
LSW theory. 

Table I The ground state properties for J' = 1, J = 0. 
Eg is the ground state energy per six spins. 





LSW 


MSW 


ED 


M A 


0.3476 


0.3955 


0.3961 


M B 


0.1951 


0.2910 


0.2922 


E g 


-2.873 


-2.910 


-2.908 



For the chain system, there is only one mean field in the 
MSW theory. It should be mentioned that neither /i a nor 
fib is zero even at zero temperature, fi a = — 0.1716J' and 
fib = 0.4001J'. It means that the solution obtained by 
the MSW theory is not a simple mean-field solution. Note 
that Eg calculated by the MSW theory is slightly lower 
than the exact E g , which may be due to the lack of the 
variational principle. For the square Heisenberg model, 
it is known that the chemical potential tends to zero at 
zero temperature [ p"4| . Thus in this case the solution by 
the MSW theory at T = is nothing but the mean-field 
solution. 

Now we proceed to the magnetic structure for finite 
J'/ J. In the limit of J' 3> J, we may use a perturbation 
theory with respect to J. It gives us SE g — 1.22OJcos0 
per unit cell, where 9 is the angle between the moments 
on neighboring ./'-chains. Thus it is seen that the quan- 
tum fluctuations favor the state 6 — tt, i.e., the stripe 
phase. For the stripe phase, the loss of the classical en- 
ergy is ~ J 2 / J' while the gain by the quantum fluctu- 
ations is ~ J. Therefore the stripe phase is expected 
for J'/J >1. The present 1/4-depleted square lattice 



2 



5=1/2 Heisenberg model may be one of the examples 
where the quantum fluctuations favor a collinear phase 

& 

Next, we turn to the other limit J' <C J. Because 
the lattice structure is bipartite in this case, the Neel 
order is expected. By use of the MSW theory, E g — 
-3.038 J per unit call, and M A = 0.2337 and Mb = 
0.2493, respectively. (fi a = 2.112 x 10~ 2 J and fi b = 
-4.090 x 10- 2 J.) 

Finally, we study spin structure for general J'/J and 
determine the phase diagram. For this purpose, we have 
to analyze the spin structure which is not stable in the 
classical sense. Thus the LSW theory is not suitable be- 
cause it produces unphysical imaginary magnon spectra, 
reflecting the classical instability. On the other hand, 
the MSW theory is applicable for such an analysis be- 
cause the self-consistent treatment of the quantum fluc- 
tuations may lead to stabilization of the spin structure 
which is unstable classically. In the MSW calculations for 
the stripe or the Neel phase, the three mean fields and 
the two chemical potentials are determined by the self- 
consistent equation and a unique solution is obtained. 
More details of the present treatments will be published 
elsewhere. 

The results of the present study are summarized in Fig. 
4 and Fig. 5 : 

(i) Stripe Phase : This structure is always unstable 
in the classical sense. Nonetheless, the quantum fluc- 
tuations stabilize this magnetic state. By use of the 
MSW theory, the solution for the ordered phase is found 
to be stable for J'/J > 0.6932. (In case J = J', 
fx a = — 0.1134J and = 0.2607J.) The recent exper- 
iments by neutron diffraction Q report that the stripe 
order is realized in CaVsO?. Note that the quantiza- 
tion axis is along the x-direction experimentally gj. In 
the present theory direction of the ordered moments is 
arbitrary because spin-orbit couplings are neglected. By 
comparing the results of theory and experiment, it is seen 
that the coupling constants J' ~ J are consistent with 
CaVsOy. Thus, we conclude that the classically unsta- 
ble (not a metastable) magnetic structure is stabilized by 
the quantum fluctuations. This is the unique feature of 
CaV 3 7 . 

(ii) Neel Phase : By use of the MSW theory, the so- 
lution for the ordered phase is obtained for < J'/J < 
0.4437. Note that the quantum fluctuations stabilize the 
Neel order beyond the classical critical point, J 1 / J = 1 /4. 
The order-disorder transition at J'/ J = 0.4437 is of the 
second order. 

(hi) Disordered Phase : The solution for the disordered 
phase is found for 0.4437 < J'/J < 0.6932. The mini- 
mum of the spin gap, which is shown in Fig. 5, is always 
at k = (0,0). The transition between the disordered 
phase and the stripe phase is of the first-order. 

(iv) Fan Phase : The LSW theory predicts that the 
energy of this phase is always higher than the others 
except for the very narrow range around ~ 0.7. (i.e., 
0.6769 < J'/J < 0.7346.) However, more importantly, 



it should be noted that the fan phase cannot be stabi- 
lized by quantum fluctuations because its spin reduction 
is always divergent. The energy for this phase obtained 
by the LSW theory may be unreliable because of lack of 
self-consistency. 

In conclusion we have shown that the magnetic struc- 
ture of CaVaOy, the stripe phase, is not stable in the 
classical theory and that only by taking quantum fluc- 
tuations into account stability of the stripe phase is ex- 
plained. For this purpose we have extended the MSW 
theory to the system where the reduction of the magnetic 
moments depends on sites. Finally we would like to em- 
phasize that the regular depletion of the square lattice 
spin-1/2 Heisenberg model realized in CaV n 02 n +i leads 
to novel phenomena, the plaquette RVB phase for n = 4, 
and the classically unstable stripe phase for n = 3. It 
is an interesting future problem to elucidate the ground 
state of the third system, n = 2, of this series. 
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FIG. 1. A model for CaVaO^ of the spin-1/2 Heisenberg 
model with the nearest neighbor (solid lines) and the next 
nearest neighbor (broken lines) exchange interactions. The 
dot dashed lines show the six-spin unit cell of the 1/4 depleted 
square lattice. 
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FIG. 3. The chain Heisenberg model constructed only by 
the next-nearest-neighbor coupling J'. 
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(a) Neel Order 
J E c i = -2 J + 2 J' 

(b) Stripe Order 
J E c i = -2 J' 



(c) Fan Structure 
( J7J > 1/4 ) 



Eel 



-2 J' 



J 2 /4J' 



FIG. 2. Typical classical magnetic structures. The fan 
phase is the ground state for J'/ J > 1/4 and the Neel phase 
is the ground state for smaller frustrations. In the classical 
treatment, the stripe phase is a saddle point and thus always 
unstable. 
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FIG. 4. Comparison between energies of different states. 
In this figure, 'Stripe-cl', 'Neel-cl' and 'Fan-cP represent the 
classical energies of the corresponding magnetic orderings. 
'Stripe' and 'Neel' represent the energy calculated by the 
MSW theory, respectively. 'Fan' represents the energy cal- 
culated by the LSW theory, whose spin reduction however is 
always divergent. Energy of the disordered state is shown by 
the thick dashed line. 
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FIG. 5. Phase diagram obtained by the MSW theory, con- 
sidering the Neel, stripe and the disordered phases. Solid 
(broken) lines represent the sublattice magnetization for 
A(B)-sites, respectively. For the disordered phase the magni- 
tude of the spin gap is shown by the dot-dashed line. The ex- 
perimental data by the neutron diffraction [2] are also shown, 
where the form factor of d xz (or d yz ) is assumed. 
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